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Abstract - The dynamics of a mean-field model 
of three coupled condensates is studied for ini- 
tial states characterized by twin condensates. 
Condensates occupy an open three-well chain, 
where the central potential-well depth is an in- 
dependent parameter. Despite its diversity from 
the closed-chain symmetric model, such model is 
shown to have an integrable regime involving a 
predictable evolution for twin-condensate initial 
states. After establishing the dependence of the 
phase-space topology on the model parameters, 
the latter are shown to allow the control of vari- 
ous macroscopic effects such as rapid population 
inversions, pulsed phase changes, self-trapping 
and long-duration periodic oscillations. Finally, 
based on quantizing the system in terms of suit- 
able "microscopic" canonical variables, the sta- 
bility conditions of the integrable regime are re- 
lated to the su(l,l) structure of the Hamiltonian 
spectrum. 
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1. INTRODUCTION 

The increasing skill in realizing experimental 
set-up's where Bose-Einstein condensates (BEC) 
are distributed in potential wells with complex 
space structures |2] has prompted, over the 
last few years, the study of models such as the 
two- well system (dimer) HJ [H], the three- 
well system (TWS) 13 |H1 E], and their many- 
well generalization such as the linear chain ^U] 
of coupled BECs. An aspect that has raised par- 
ticular interest is the recognition of characteris- 
tic behaviors and macroscopic phenomena that 
distinguish the dynamics of coupled wells. The 
latter consists of interwell boson exchanges en- 
acted by the tunnel effect that determine, in 
each well, macroscopic changes of the condensate 
phase and population. 

The study of the TWS on a closed chain (three 
wells mutually coupled) has revealed, within a 
dynamics dominated by chaos the presence 
of an integrable subregime |H1 El corresponding 
to (a phase-space submanifold of) states in which 
two of the three wells exhibit identical conden- 
sate phases and boson numbers per well. Such 
twin-condensate (TC) configurations are pre- 
served dynamically and make the TWS equiv- 
alent, in practice, to a two-well system. Partic- 
ularly, the dimer integrability [HI is inherited by 
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TWS when this is prepared in an initial state 
with TCs. 

The TWS integrable subdynamics raises in- 
terest for various reasons. First, achieving a full 
knowledge of the TWS phenomenology with a 
regular dynamical character (this extends the 
analysis of Ref. is necessary to design sig- 

nificant experiments. Second, since perturbing 
dimerlike states leads to prime chaotic behaviors, 
the modalities governing their onset (namely 
the preparation of suitable initial states priming 
chaos) are a central aspect. This is important, 
not only for the TWS, but also for the new class 
of systems (coupled-BEC chains) whose control, 
at the experimental level, is rapidly increasing. 
Finally, in view of the last observation, the recog- 
nition of a class of states having regular time be- 
haviors (near those showing chaos) is certainly 
useful to study the manifestation of classical in- 
stabilities at the quantum level. 

In this paper we consider the open TWS, 
namely a three-well system in which two decou- 
pled (lateral) wells interact with a third (central) 
one. In particular, we show that, as in the closed- 
chain case, the open-TWS dynamics is equipped 
with an integrable subregime. To do this, we find 
new coordinates that both account for the dy- 
namical constraint embodying the TC form and 
preserve the canonical structure in the TC sub- 
dynamics without resorting to Dirac's procedure 
for constrained Hamiltonians. As to the possibil- 
ity to influence the system dynamical behaviour 
(and thus to play a relevant role in designing pos- 
sible controlled experiments), we point out that 
in addition to the tunneling parameter (denoted 
by T) the open TWS model exhibits a further 
adjustable parameter w that distinguishes the 
central potential-well depth from the two (equal) 
lateral- well depths. 

The present analysis prosecutes the work both 
of Ref. 9 concerning the self-trapping onset in 
the dimerlike regime of the symmetric {w = 0) 
closed TWS, and of Ref. In the latter 

the analysis of the fixed points and the origin 
of chaotic behaviors in the open TWS were em- 



bodied, for arbitrary T and w, in two dynamical 
phase diagrams (relevant to the dimerlike and 
nondimerlike classes of fixed points) that supply 
an exhaustive, operationally valuable, account of 
the system stability. 

In the sequel, we concentrate on the dynam- 
ics of TC configurations and, more specifically, 
on the features "in large" of orbit bundles that 
characterize the topology of the reduced phase 
space in the dimerlike sub-regime. We investi- 
gate both the conditions under which the sys- 
tem evolution exhibits either oscillatory or bal- 
listic motion, and the possibility to control the 
self-trapping and other effects through the pa- 
rameters T and w. Also, we show that the inde- 
pendence of parameter w, which introduces fur- 
ther asymmetry in the model-Hamiltonian struc- 
ture but improves the possibility of fitting ex- 
perimental conditions, does not affect the (par- 
tial) integrability of TWS and enriches notice- 
ably the dynamical scenario of the system. Fi- 
nally, we evaluate the stability of TC dynamics 
under quantum perturbations by partially recov- 
ering the quantized form of the TWS Hamilto- 
nian. To this end the dynamical degrees of free- 
dom that, within the new canonical coordinates 
mentioned above, account for the deviation from 
dimeric regime, naturally undergo the quantiza- 
tion when they have a nonmacroscopic character. 
The stability of dimerlike configurations will be 
proven by relying on the su(l,l) structure exhib- 
ited by the model Hamiltonian. 

2. TRIMER-DYNAMICS CANONICAL FORM 

The asymmetric TWS we consider here is de- 
scribed by the second-quantized Hamiltonian 

3 

H = U^^ nf - v{ni + n^) - V2n2 

i=l 

T 

- — {a\a2 + ajas + h.c)j , (1) 

which one derives from the many-body quantum 
theory of BECs through a three-mode expansion 
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of the condensate field operator Parame- 
ters U, T, V, account for the interatomic scat- 
tering, the tunnehng amphtude, and the sided- 
well potential depth, respectively; the central- 
well depth has, in general, an independent value 
V2 = V + w ^ V. Another element of diversity 
from the symmetric-TWS case is due to the ab- 
sence of interwell coupling between the first and 
the third well. The operators ni = a\ai count 
the bosons in the rth well (A^ = Sjnj), while the 
destruction (creation) operators (a|) obey the 
canonical commutators [aj,aj] = da. Preceding 
studies of the TWS dynamics have been focused 
on a case where the asymmetric character is- 
sues from nonconstant tunneling parameters Tij 
such that Ti2 ^ T13 = T23 = T. Classically 
{a\ = a|, OjoJ = a\ai), the asymmetric TWS 
has revealed |i7j| the presence of homoclinic chaos, 
while, quantally, the survival of breather con- 
figurations has been investigated on the trimer 
viewed as the smallest possible closed chain. The 
Heisenberg equations related to H for the boson 
operators Oj, a\ give, within the random-phase 
approximation, the equations (s = 1,3) 



ihzs = {2U\zs\ - v)zs - ^Z2 



ihz2 = [2f7|22p -{w + V)]Z2 - ¥(^3 + Zi) 



for the expectation values Zi = (ai), z* = {a\), 
of the three wells. These entail Ej|zjp as a con- 
served quantity replacing the (conserved) total 
boson number N {[N,H] = 0). Eqs. @ can be 
also obtained from 

W(Z,Z*)^t/f]|z,f -t;(|zi|2 + |z3|2) 

-(t; + ^11)122^ - '^{zlz^i + zlzi + c.c.) 

by using the standard canonical Poisson brackets 
{z1,,Zj} = i5kj/h. The latter and 7i{Z,Z*) can 
be also derived variationally within the coherent- 
state approach developed in Refs. 1^ . 

Sub-dynamics of Twin-Condensate config- 
urations. Based on Eqs. ((21), imposing the con- 



straint zi = Z3 leads to identify an integrable 
subregime of TWS dynamics which consists of 
the orbit bundles lying on the sub-manifold V := 
{(2:1,2:2,2:3) : zi = 23,2:2 S C} in the original 
phase space A4. This can be proven by introduc- 
ing the new canonical variables ^ = (21 — 23) /V2, 
z = [zi + 23)/\/2 such that {d] = ijfi = 
{2*, 2}, and {2*,e} = {e,Z2} = {z*2,z} = 0. 
One finds that 



H 



u_ 
2 



{\Z\' + \i\'Y + {zC + Z*iY + 2\Z2\ 



vM - 



'^i\zi\ 



W 22 



\Z2? + \Z\'^ 



V2 



iz2Z*+z*z) (3) 



giving, in turn. 



ihi = U{\z\^ + \^\^)z-vz+U{zC+z*0^-^Z2 



ihi2 = [2U\z2\ 



{v + w)]z2- -^z, (4) 



ihi = C/(|2|2 + |C|2)^ - < + UizC + Z*i)z . 

Setting ^ = (that is 21 = 23) eliminates the 
third equation, while the remaining pair of cou- 
pled equations, with the substitution 2 = \f2z\^ 



(2) 



ihix = (2f/|2i|2 

ini2 = [2[/|22|2 ■ 



V Zi 



T 

2^2 



(5) 



{v + w)]z2 - Tzi 



coincide with Eqs. under the restriction 21 = 
2:3 to the sub- manifold V. The two costants of 
motion corresponding to the energy 



U [ ^W + \Z2\ 



v\z\ 



-iv + w)\z2\'^ - ^{Z2Z* 



Z2Z) 



and the total boson number N = 2ni + 71-2 



kl + 1-221 i^i 



make the dynamics issu- 



ing from Eqs. © integrable. Notice that, while 
Hq still generates the correct equations for z and 
22 through 2 = {2, Hq}, 22 = {22, -ffo}) the same 
Hamiltonian written in terms of 21 , Z2 no longer 
provides Eqs. ©. The useful role of the vari- 
ables 2, 22, ^ is thus to supply a new canonical 
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scheme where the restriction to V can be enacted 
directly on H bypassing Dirac's scheme for con- 
strained systems. The resulting Hamiltonian Hq 
can be further cast in the pendulum-like version 



Ha 



3U 



fu 



cos 26 . 



in which the motion constant is used explic- 
itly, Co = 3UN^/8 - N{v + w/2), = UN{1 - 
2v), V = w/UN, and 6 := {(1)2 - (p)/2 [defined 
through z = |z|exp(i0), Z2 = |z2| exp(z02)] and 
D := |z|2-|z2p fuimithe brackets {6*, D} = l/h. 
The ensuing representation on the plane {9, D) 
is the simplest possible way to describe the space 
V and its structure changes when T, U, w are 
varied. Such a 9 — D picture involves two equa- 
tions 

TD cos 29 



he = -UD - 4^ + ^ ^ 

4 4 V2^/W 



J)2 ' 



nt) = -V2tVn^ - sin 29 . 



(6) 



(7) 



only. The other pair of canonical variables ip = 
{(p2 + (/')/2, N = |zp + |z2p, in fact, do not 
participate in the dynamics: while = 0, the 
angle ip plays an auxiliary role since ip can be 
proven to be completely determined by the evo- 
lution of 9 and D which, on the contrary, have 
V'-independent equations. 

3. PHASE-SPACE STRUCTURE 
AND PARAMETER DEPENDENCE 

The TWS integrable subdynamics is reducible 
to a one-dimensional potential problem with 
an unique (phase-independent) E'A'^-dependent 
equation = f{D;E,N) for D by implement- 
ing standard quadrature procedures [HI E] • The 
latter lead to solve explicitly the motion equa- 
tions in terms of elliptic functions. Nevertheless, 
the analytic complexity of the one-dimensional 
potential makes rather difficult to establish the 
dependence on r and of the V structure as well 
as of the dynamics therein represented. Then, 
since we wish to recognize the features that dis- 
tinguish the dynamical behaviors of TC configu- 
rations [rather than the exact analytic solutions 



of Eqs. (ini), 0], we simply approach the prob- 
lem by performing in parallel the study of the 
fixed points of Eqs. ©, and the numeri- 
cal reconstruction of the phase-space portraits. 
This furnishes in a quite direct way information 
on both the topology changes in V when r and 
v are varied, and the features "in large" of the 
dimerlike bundles of trajectories represented on 
the 9 — D plane. In particular, the conditions 
that cause variuos macroscopic effects in the dy- 
namical behavior can be evidenced effectively. 

The equation pair that determine the fixed 
points, = [WD-U+2V2TD cos{29) /r{D)]/4 , 
= V2Tr{D) sm{29) with r{D) = VN^ - D'^, 
is stemmed from equations ©, © when setting 
9 = = D. 

These supply maxima, saddles and minima 
that constitute the fixed-points set in the re- 
duced phase space V. The minimum is situated 
in 6* = 0, L> / (L> = only for u = 1/2), while 
the other extremal points stay on the straight 
line 9 = tt/2. Under such conditions, the val- 
ues of D at fixed-points are obtained from the 
nonlinear equation 

= ^X-1^+^^I^ (8) 
4 4 V2VT^^ 

where a = +1 (cr = —1) in the case of the min- 
imum (maxima/saddles), and the adimensional 
quantities r := T/UN, X = D/N have been in- 
troduced in addition to := w/UN. The num- 
ber of fixed point function of i/, r is easily 
deduced from eq. ©. In the parameter space 
(i/, r), the critical lines 



2^ ' 2 



2/3' 



-1 3/2 



(r = 1,2) are recognized to play a crucial role. 
For 1^1 (r) < V < i'2{'t) = ^^i (''")) 2^2 ("?")) three 
(two) solutions are furnished by equation © 
with (T = — 1, in addition to the permanent single 
solution of the case a = +1. The lobe 

C = {{u,t) : < r < 3/^8, ui{t) <v < V2{j)), 
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defined in the parameter plane (z^, r) identifies 
completely the area corresponding to the three- 
solution regime. Outside £, only one solution 
survives together with the solution of the case 
a = +l. 

4. TOPOLOGY CHANGES AND 
DYNAMICAL EFFECTS 

The recognition of the critical lines bounding 
C together with some numeric simulations of the 
reduced phase space V (see the figure pairs 1 and 
2), enable us to determine the influence of r and 

V on the topology of 7^, namely on the structure 
of orbit bundles of the TC dynamics. 

For (i/, r) G £, the fixed points resulting from 
Eqs. 0, with cr = +1, consist of two maxima 
M_, M+ and a saddle Mq such that 6* = 7r/2 
and X_ < Xq < Xj^. Such a situation ^Jj is 
described in Fig. la and Fig. lb where, with 
r = 0.5, the cases v = 0.25 and u = 0.75, respec- 
tively, are depicted. Their symmetric structure 
under the global reflection with respect to the 
axis D/N = 0, reflects the symmetry of Hamil- 
tonian i^o which combines the changes v ^ 1 — u 
and D -D. 

When r = const and u ^ (u ^ i^^), a 
coalescence effect takes place in which M_ (M_|_) 
and Mq merge in a unique point. This process 
can be driven either by the potential-depth gap 
w = ^2 — f or by the total boson number 
through the parameter u = w/UN (notice that 
U, accounting for the s-wave scattering length, 
in the present paper is considered as a fixed pa- 
rameter). The presence of the adjustable gap w, 
however, allows one to operate leaving (and 
thus r = T/UN) unchanged. Fig. 2a and Fig. 
lb illustrate the effect of crossing the C border 
by changing u from 1 to 0.75. No change of the 

V structure occurs by further decreasing v as far 
as i'i{t) < V. Fig. la and Fig. lb, in fact, have 
the same number of extremal points. 

Thanks to the presence of r, the coalescence 
can be achieved as well if, keeping u constant 
within the interval — 1 < < 2, the change r 
r' > r is such that (z/, r) G £ while {v, t') 



C. The comparison of Fig. la and Fig. 2b, 
both having v = 0.25, well illustrates the phase- 
space structure emerging from the coalescence of 
M+ with the saddle Mq as a consequence of the 
change r = 0.5 — > r' = 1. 

The orbit-bundle structure deeply reflects the 
action of z^-r changes. For {u, t) £ £, V dis- 
plays five independent bundles whose topology 
is characterized by means of three separatices 
(dashed lines in all the figures) 5+, 5_, and Sq 
(the latter shows a 7 shape). Outside C, owing 
to the coalescence (see Figs. 2a, 2b), separa- 
trix Sq disappears, while the orbit bundles re- 
duce to three . Curves S± are easily identified by 
observing that they start/terminate at D = ±N, 
6 = 7r/4, Sir/ 4, whereas Sq is based at D = NXq, 
6 = tt/2, where Xq is given by Eqs. (ISJ. Notice 
that for 1/ < 0.5 the two basins encircled by 
and S- contain the minimum m and the maxi- 
mum Af_, respectively (see Figs, la, 2b), while 
the opposite situation appears for v > 0.5 (see 
Figs, lb, 2a), where 5+ and S- encircles M_|_ and 
m, respectively. In both cases, such basins are 
filled by closed-curve bundles entailing bounded 
oscillations for 6. 

Two bundles of ballistic orbits, instead, are com- 
prised between 5+ and 5_ and separated by Sq 
when (z^, t) G C. In this case, no constraint lim- 
its 9 that covers its whole range [0, 2tt], as shown 
by Figs. 3a, and 4a. The closed arc Co, consti- 
tuing the noose of Sq, confines a further closed- 
orbit bundle that encircles M_|_ (M_) for u < 0.5 
{v > 0.5). 

When the coalescence is enacted, Co and the 
basins around M± collapse, while Sq merges to 
the ballistic-curve flows lying between S- and 5+ 
(see Figs 2a, 2b). This topology change is accom- 
panied by another, dynamically relevant, effect 
that consists of the suppression of the two mo- 
tion curves forming Sq: the closed branch Co and 
the A-shaped branch. Such orbits, both starting 
and terminating at Mq, are characterized by an 
infinite covering time (illustrated in the follow- 
ing sub-section) as required by the saddle fea- 
ture of Mq. Orbits placed in their proximity. 
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which before the coalescence display arbitrar- 
ily long percurrence times, lose such a character 
after the saddle suppression. A special macro- 
scopic trait therefore distinguishing the param- 
eter choices (i^, t) £ C and {v, t) ^ C is that 
the former choice always involves the presence 
of nonzero set of initial conditions around Mq 
whose characteristic time period can become ar- 
bitrarily large. This feature, of course, is inter- 
esting both for achieving a complete control of 
the TWS characteristic times and in view of the 
experimental observations. 

Ballistic motion and self-trapping. In Figs. 
3, where the elapsing time is UNt, we illustrate 
the evolution of both D/N (its oscillations are 
evidenced in the grey stripe) and 6 relevant to 
various orbits in 7^. In Fig. 3a, D/N and 9 are 
given for two orbits having energy values Ea, 
Eb, with T = 0.5, z/ = 0.25. These energies 
are such that Ea — Eb — 5 • W~'^Eo, where Eq 
is the energy of Sq and \Eb — Eq] ~ 10~^£'o. 
The period of the Eb orbit is greater than that 
of the Ea orbit since the former is closer to 
Sq. Both 6 and D/N exhibit an intermittent 
evolution where time changes are confined in a 
sub-interval of the characteristic period in which 
rapid population inversion effects are visible as 
well. Fig. 3b shows D/N and 6 for two orbits 
near m, and M_ and allows one to compare the 
corresponding periods with the ballistic motions 
of Fig. 3a. The evolution of and D/N is dis- 
played in Fig. 4a for an orbit almost coinciding 
with So (rapid local variations of both are visi- 
ble in correspondence to the motion around the 
So noose). Fig. 4b shows and D/N for two 
orbits near 5_ for t = 1, u = 0.25: the inner or- 
bit exhibits a pulsed population inversion where 
abrupt changes alternate with an almost linear 
decreasing of D/N. 

In consequence of the two-parameter depen- 
dence of the V topology, also the self-trapping 
effect can be induced/suppressed through r and 
v. We recall that the latter is particularly evi- 
dent IHl El when the phase space displays pairs 



of isoenergetic orbits. In this case the system is 
imposed to single out one of the available orbits 
at a given energy. This symmetry breaking effect 
entails the self-trapping in a restricted region of 
V which cannot be left unless sufficiently large 
energy changes are effected. 

In the presence of two, well separated max- 
ima M_ and M+ (Fig. la and Fig. lb), the 
self-trapping regimes exhibit oscillations of pop- 
ulations Hi = \zi\'^ such that either ni = ns ~ 0, 
n2 ^ N (orbits close to M_ with D ~ —N) or 
JT-i = "fia — n2 — (orbits close to M+ 

with D ~ N). As a result, the trajectories encir- 
cling the maximum manifest an evident, stable 
population imbalance (namely ni = » n2, 
ni = n3 « n2) during the system evolution that 
excludes the oscillations to approach intermedi- 
ate states rei = ns ~ N/4, n2 ~ N/2 {D ~ 0) 
near the minimum energy configuration m. We 
wish to notice that, when the suppression hap- 
pens and M_ (-M_|_) and Mq merge in a unique 
point, in general, the remaining maximum M_|_ 
(M_) keeps a position that is rather different 
from that of the minimum (see Fig. 2a, Fig. 2b) 
and thus preserves a regime with a marked pop- 
ulation imbalance. 

Concerning the change v ^ 1 — v (<^ D — > 
—D) it is worth noting that, while the basins of 
Mo and M± remain unaltered, the ballistic orbit 
bundles between and S- undergo an abrupt 
inversion. For v = 0.5 the linear D term in Ho 
disappears and the resulting symmetry D —D 
renders the dimeric regime of the open TWS 
equivalent to the pure symmetric dimer 
We devote a final remark to the case i/ = (equal 
well depths). Since (r, i/) £ for any r 7^ if 
I' = 0, then V portraits are quite similar to the 
one of Fig. 2b showing the ballistic phase evolu- 
tion and pulsed population inversions (see Fig. 
4b). It is worth noting that the evident asym- 
metric character of the TWS dynamics, well dis- 
played by the V portraits, ensues not only from 
the well-depth diversity but also from the fact 
the pair of TC is described within the integrable 
regime by an unique collective variable z. 
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Figure 1: (a) Left figure is obtained by setting r = 0.5, v = 0.25; (b) right figure corresponds to 
T = 0.5, V = 0.75. In both figures, the minimum m, the saddle Mq and the maxima M± are visible 
as well as separatrices Sq, S± represented by different dashing styles. 
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Figure 2: (a) Left figure illustrates V ior t = 0.5, v = 1: increasing u from v = 0.75 causes the 
coalescence of Mq and M_ displayed in Fig. lb. (b) Right figure corresponds to t = 1, u = 0.25: 
increasing r from the value r = 0.5 causes the coalescence of Mq and M_|_ displayed in Fig. la. 
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Figure 3: Evolutions of D/N G [—1, 1] (in the grey stripes) and 9 for some orbits in V (see the 
sub-panels), with r = 0.5, v = 0.25; UNt is the rescaled time, (a) Left figure: 9 and D/N for two 
orbits very close to the A branch of Sq. The orbit closer to Sq has a larger period, (b) Right figure: 
periods of left-figure orbits can be compared with those of two orbits near M_ and m. Refer to 
the text for further comments. 




Figure 4: Evolutions of D/N G [—1,1] (evidenced in the grey stripes) and 9 for various orbits in 
V (see the sub-panels). UNt is the rescaled time (a) Left figure corresponds to an orbit very near 
»So, with r = 0.5, u = 0.25. (b) Right figure shows two orbits adjacent to S- with t = 1, v = 0.25. 
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5. QUANTUM APPROACH TO 
DIMERLIKE DYNAMICS 

The trimer dynamics naturally undergoes a 
quantum formulation both when it is investi- 
gated in the low-energy regime and if small par- 
ticle numbers per well do not allow Bogolubov's 
approximation. Classically (namely in the mean- 
field picture described in section 2), the ground- 
state configuration is fully determined and rep- 
resents a special case within the dimeric regime. 
In the alternative canonical scheme discussed 
above, quantum Hamiltonian takes the form 
(ric := c^c, c = a,b) 



H = U 



{ua + ribf + \{h^a + a)hf 



first order terms in ^, Under such assump- 
tion the third equation, where is now neg- 
ligible, is independent in that the time evolution 
of z-dependent parameters is determined by the 
first two equations. Also, its form implies that 
1^1 << |i|, \z2\ thus making ^ a slowly varying 
variable. This scenario suggests that the quan- 
tum problem for the ^ mode can be described by 
the Hamiltonian 

Hi, = ^{2\zfnh + nl) + ^{hU + z*hf - vnt , 

where is negligible with respect to 2\z\'^nb. 
The spectrum of Hb is easily calculated by recog- 
nizing that its spectrum generating algebra ^H] 
is su(l,l). The latter is generated by 



T 

-v{na + Ub) - wn2 - -j= {a2a} + aja) i i^s = (2nb l)/4, = 



where the quantum counterparts of new vari- 
ables z are 

6 = (oi - as 



2, 



(ai +a3)/V2, 



obeying the commutators [6,6^] = 1 = [a, a"^]. In 
the spirit of Bogolubov's approximation (which 
is valid when a bosonic mode c is strongly oc- 
cupied so that ~ c* and c^c ~ |cp with 
[c, c^] ~ ), we observe that for states close to 
the dimeric regime both a and 02 are macro- 
scopic, while 6, 6^ can be seen as microscopic de- 
grees of freedom describing the deviation from 
the pure dimeric states. In particular, the clas- 
sical condition ^ ~ should correspond quan- 
tally to have the expectation value (6) ~ 0, 
whereas -classically, the population of the 
^ mode- should correspond to the operator rib- 
Only quantum states involving small quantum 
numbers for Ub are expected to be excited in the 
system dynamics. 

In the light of this, Eqs. (jlj) can be employed 
to approach in the appropriate way the dynamics 
in which b, and 6^ are quantum objects, whereas 
a and 02 can be identified with z and Z2, re- 
spectively, owing to their macroscopic charac- 
ter. The first two equations are in fact decou- 
pled from the third one if one considers only the 



with commutators [K^,K±] = ^K±, [K-,Kj^] = 
2K3. Upon noticing that (J}^z + z*bf = z^{b^f+ 
(z*)^6^-|- |2;p(2nf, -|- 1), one gets Hb in the new 
form 



Hb = 2(2C/|z|2 - v)K3 - y 



2 1^1 



+U[z'K+ + {z*yK^] + - . 

By using the unitary transformation U3 = 
ex.p[i(j)K3], (notice that U3K±ul = e^''l'K±), 
where (/> is given by z = \z\e^'^, one finds that 
Hb = U3HM with 



H^ = 2{2U\z\^ - v)K3 + 2U\z\'^Ki 



z + 



2 ' 



where Ki = (K^ + K_ )/2 is one of the two non- 
compact generators of the pseudo-angular mo- 
mentum standardly defined by means of K± = 
Ki lb 1X2- The further transformation U2 = 
exp[— iaii'2] such that 

U2K3UI = di{a)K3 + sh(a)Ki , 

allows H^ to be expressed as 



H, = U2[RK3 + -'^\z\'' + l]Ui 
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= Rch{a)K3 + Rs\v{a)Ki - ^IzP + - 

^ 2 

provided tha = U\z\'^ / {2U\z\'^ - v) and R^ = 
(3C/|2;p — v){U\z\'^ — v), and either |zp > v/U or 
jzp < v/'iU is fulfilled. Whenever one of the lat- 
ter conditions is violated {v/3U < [zp < v/U) 
the transformation U2 must be related to the 
noncompact element Ki (instead of K^) in or- 
der to reproduce H^:. The ensuing effect is dra- 
matic: While K3 has a discrete, positive spec- 
trum bounded from below, 

spect{K3) = {2n + l)/A, n = 0,1,2,..., 

Ki has a continuous, unbounded spectrum. 
Since we have used unitary transformations such 
spectra are those characterizing Hi, itself. There- 
fore, only in the case when the initial conditions 
on z ensure the presence of a discrete spectrum 



(2n + 1) 



En{\z\^) = ^mA''-v){U\z\^-v) 



u 



the ground-state turns out to be stable in that 
an energy gap separates it from the first excited 
state. This makes stable the dimeric regime. A 
final check is quite natural at this point. Making 
use of the first of eqs. Q in the integrable sub- 
dynamics = 0) and recalling that z = \f2z\, 
the discreteness condition becomes 



4C/|zir + 



TZ2 \ TZ2 



2 ziJ 2 zi 



> 



(9) 



Remarkably, the conditions on the Hamiltonian 
parameters r and i' enforcing condition ® are 
basically the same ensuring that the classical 
fixed points are stable. That is to say, the dis- 
creteness condition Q selects the stable regions 
(maxima and minima) of the dimeric stability 
diagram appearing in ref . ■ 

We briefly recall that, owing to a global phase 
symmetry, each dimeric fixed point can be iden- 
tified by a pair of real numbers, xi and X2, such 
that zi = Z3 = \/iV xi and Z2 = X2. It 



proves then useful to introduce the real param- 
eter Q = X2/{V2xi) such that xi = cos(i?)/\/2, 
X2 = sin(i?), where 'd = arctana. This way each 
fixed point is uniquely related to a single real pa- 
rameter, either or a = tan It turns out that, 
for a given choice f,^ of the parameters, the rele- 
vant t?'s can be found by intersecting the straight 
line T = f and the curve r = max[0. To {{))], 
where r^(t?) = [u + I + a^{v - 2)]/(a^ - 1). 
If ■!? > the relevant fixed point is a mini- 
mum. Conversely, if < 0, it is either a max- 
imum or an unstable saddle depending on 
whether f > f[{}) or < f < f{'&), where 
f(i?) =max[Ti(^?),r2(^?)], 



ri{e) 



-2^/2 



/3(l + /?2) 



, T2{e) 



-6V2P^ 

(l + /?2)3' 



(10) 



with P = tanO. Now, recalling that T = UNt, 
equation ® finally reduces to 

[2r tan 61/(1 -htan^ 61) -hT2tan2 6l/\/2] > (11) 

which selects the minima and, almost exactly, 
the maxima regions defined by the study of the 
stability character of the classical fixed points 
above recalled. Indeed the points in the minima 
region, {0 < 'd < tt/2,t > 0}, trivially meet con- 
dition (|llj) . The same is true for most of the 
maxima region, {—tt/2 < ?? < 0, r > f{i))}. In- 
deed, if 1? < 0, condition (fTTI) selects the region 
above ti(t?), which is exactly the borderline be- 
tween maxima and saddles when ti{9) > T2{9). 
In the complementary case, ti{6) < T2{6), a 
small portion of the saddles region is predicted 
to give rise to a discrete rather than continuous 
spectrum. This slight disagreement can be ten- 
tatively explained by observing that in the lat- 
ter case "d (is likely to be) can be quite close to 
— 7r/2, which ultimately means \zi\ ~ 0. This, 
however, is in contrast with our hypothesis that 
zi is macroscopic. 

5. CONCLUSIONS 

Based on a variational technique previously 
developed [HI we have studied the integrable 
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sub-regime of the TWS dynamics in terms of the 
expectation value Zi (z*) relevant to operators 
Oj (a|). Our approach is equivalent to the stan- 
dard mean-field picture. Upon defining suitable 
collective canonical modes, in section 2 we have 
presented an alternative description of the dy- 
namics which allows to project the TWS dynam- 
ics on the reduced {zi = z^) space V containing 
the integrable orbits, and avoids the nontrivial 
Dirac's method for constrained hamiltonian sys- 
tems. The new picture thus obtained, in which 
TCs are always present, turns out to be equiv- 
alent to an asymmetric dimerlike model whose 
integrability is well known. 

After deriving the fixed points of Hq and the 
conditions determining the coalescence of fixed- 
point pairs depending on r and v (section 3) , we 
have studied, in section 4, the TWS integrable 
sub-dynamics. The changes of the V topology 
and, more specifically, the changes of the struc- 
ture of orbit bundles filling V have been consid- 
ered under the variations of r and z^, and the en- 
suing effects on dynamics have been dicussed. In 
particular, the ballistic orbit bundles have been 
found to have a periodic evolution whose time 
scale increases when the orbit approaches the 
separatrix Sq. The possibility of having arbitrar- 
ily large time periods for a nonzero set of orbits 
disapperas as soon as a coalescence process de- 
stroys Sq. Various situations in which the phase 
undergoes either a pulsed evolution or a stepwise 
variation while D/N shows rapid population in- 
versions have been discussed. These effects share 
the common feature of being macroscopic, offer 
new information on the special TC regime of the 
three-well system, and thus are interesting in or- 
der to design new experiments on coupled BECs. 
The analysis we have performed on the V topol- 
ogy extends the study of the two-well dynam- 
ics 15, and recognizes the presence of the two- 
well phenomenology within the TWS dynamics. 
We notice that similar phenomena have been ob- 
served in the coherent-mode dynamics of weakly 
excited condensates |22] , the motion equations of 
which exhibit in fact a rather similar nonlinear 



form. 

We emphasize the fact that, beyond the rich 
scenario of dynamical behaviors, our analysis 
clearly displays both the way to prime macro- 
scopic effects and the circumstances in which 
these are expected to occur by changing in a 
controlled way the relevant parameters N, r and 
ly. In section 5 we have presented a preliminary 
study of the purely quantum TWS. A nice result 
which confirms the utility of the new canonical 
picture adopted in section 2 is that the algebraic 
analysis of the microscopic degrees of freedom 
capable of perturbing the integrable subregime 
permits to evidence the stable nature of TC con- 
figurations. These, in a sense, appear to be pro- 
tected by the energy gap ensuing from the dis- 
crete form of the energy spectrum in the compact 
sector of the hamiltonian algebra su(l,l). No- 
ticeably, the noncompact sector giving a contin- 
uous spectrum turns out to characterize the sys- 
tem in correspondence to unstable TC states. A 
systematic quantum study of the TWS along the 
lines of Ref. which is presently in progress, 
will be presented elsewhere. 
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